Theoretically predicted picosecond optical switching of spin chirality in multiferroics 
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We show theoretically with an accurate spin Hamiltonian describing the multiferroic Mn per- 
ovskites that the application of the picosecond optical pulse with a terahertz frequency can switch 
the spin chirality through intensely exciting the electromagnons. There are four states with differ- 
ent spin chiralities, i.e. clockwise and counterclockwise ab/bc-plsme spin spirals, and by tuning the 
strength, shape and length of the pulse, the switching among these states can be controlled at will. 
Dynamical pattern formation during the switching is also discussed. 

PACS numbers: 75.80.+q, 75.85.+t, 77.80.Fm, 75.10.Hk 



Chirality, i.e., the right- and left-handedness of struc- 
ture, is one of the fundamental concepts penetrating 
through the whole of science. In solids, electron spins 
sometimes form a chiral order, which offers an opportu- 
nity to manipulate the chirality by external parameters. 
This issue is of vital importance in spintronics, which 
aims at the electric control of spins by an electric current 
or electric field (E). Multiferroics provides us an ideal 
system for this purpose, in which the spin chirality is di- 
rectly related to the electric polarization [H, 0] . In the 
spin-current model [H 0] , two mutually canted spins S% 
and Sj generate the polarization as 




Pij CX &ij X (Si X Sj), 



(1) 



with e-ij being the vector connecting ith and jth sites. 
Here the vector product Si x Sj is the spin chiral- 
ity, which characterizes a direction of the spin rotation. 
In multiferroic Mn perovskites such as TbMnC>3 and 
DyMnOs, the Mn spins are rotating within the be plane 
to form a cycloid (fee-plane spiral) propagating along the 
+b direction (Pbnm setting) as shown in Fig.QJa) p, 0, 
and the clockwise (CW) one with (Si x Sj)\\—a in- 
duces P||+c, while the counterclockwise (CCW) one with 
(Si x Sj)||+a induces P\\— c. 

In this Letter, we theoretically propose a picosecond 
optical switching of the spin chirality in TbMn03 as 
a representative material. However, the mechanism is 
rather general and would be ubiquitous in spiral magnets. 
It is found that the oscillating E component of the opti- 
cal pulse activates collective rotations of the spin-spiral 
planes via magnetoelectric (ME) coupling, and their in- 
ertial motions result in the chirality reversal or flop. By 
tuning the strength, shape and length of the pulse, the 
spin chirality is shown to be controlled at will. 

When the electric polarization is driven by the spin or- 
der, it is naturally expected that the collective excitation 
of spins (magnon) has an infrared activity 0, HI- Indeed, 
strong optical absorptions were experimentally observed 
in RMn0 3 (R=Tb, Dy, Eui.J^ etc) at THz frequen- 
cies, and they were ascribed to magnons activated by 
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FIG. 1: (Color online) (a) Spin configuration, spin chirality 
Si x Si+b and ferroelectric polarization P in the clockwise 
6c-plane spin spiral (upper panel) and those in the counter- 
clockwise one (lower panel), (b) Modulation of the in-plane 
ferromagnetic exchanges under U||±a. Upper (lower) signs 
in front of AJ a b (>0) correspond to the modulations under 
-E||+a (.E||— a), (c) Calculated electromagnon spectrum. In- 
set shows the experimental spectrum for TbMnOa 



the E component of the light, i.e., electromagnons 0,0- 
14[ . In the early stage, the corresponding magnon modes 
were interpreted as rotation of the spin-spiral plane with 
oscillating in Eq. (1) @. However, this interpreta- 
tion contradicts the experimental observation about the 
selection rule in terms of the light polarization 
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Afterwards, it turned out that the polarization due 
to the conventional magnetostriction, i.e., 



Pij = *ij(Si ■ Sj) 



(2) 



is 

i?MnO 



relevant to 



the infrared-activity of magnons in 
Here the vector TZij is nonzero because 
of the orthorhombic lattice distortion without inversion 
symmetry at the center of the Mn-O-Mn bond. The puz- 
zling electromagnon optical spectrum with two specific 
peaks was successfully explained by this mechanism [l7| ■ 
Under this circumstance, the photo-induced phenom- 
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ena become a challenging issue. Since p^ in Eq. (2) does 
not require the spin-orbit interaction, its magnitude is 
much larger than that of in Eq. (1), which enables the 
intense and fast optical excitation of magnons. This of- 
fers a unique opportunity to study the nonlinear dynam- 
ics of the spin system. In addition, the light can locally 
activate or modify the spin structure with a squeezed 
light spot in contrast to the magnetic field. 

We are now ready to attack such phenomena in 
i?Mn03 theoretically ahead of experiments for the fol- 
lowing reasons. First, we know that the optical pulse ac- 
tivates mostly the spins only via the ME coupling, which 
allows us to neglect electronic excitations at much higher 
energies (>1.5 eV). Second, we have an accurate spin 
Hamiltonian, which describes competitions among vari- 
ous phases in i?MnC>3 [HI , so that the optical switchings 
among them and dynamics after the light irradiation can 
be simulated in a reliable way. 

We employ a classical Heisenberg model on a cubic 
lattice, which contains not only the frustrating spin ex- 
changes ("Hex) but also the single-ion spin anisotropy 
("H^ a and Hf ia ), Dzyaloshinskii-Moriya (DM) interaction 
(%dm) and biquadratic interaction (Hbiq) as, 

T~L = Hex + T~L^ a + %si a + Hdm + %t>iq 



dition. For the ME c oupling 
Pij given in Eq. (2) 



IT] 



= J2 JijSi-Sj + D^Sh 

<i,j> i 

i 

ab 

+ d * ■ & x Sj) - B hiq (Si ■ Stf. (3) 



<hj> 



<hj> 



The frustration between ferromagnetic (FM) exchange 
J a b and antifcrromagnctic (AFM) exchange Jb results 
in the in-plane spiral spin orders, while the interplanc 
AFM exchange J c causes their staggered stacking. The 
DM vectors dij are expressed using five DM parameters, 
a a b, Pab, lab, a c , and /3 C , as given in Ref. UJ- Crucial 
roles of the biquadratic interaction in i?MnC>3 were un- 
covered in recent theoretical studies 17, 2(J. For more 
detail of the model, see Ref. [17J. We adopt the fol- 
lowing parameters: J a b=— 0.74, Jb=0.64, J c =1.0, (a a b, 
Pa,b, 7ab)=(0.1, 0.1, 0.14), (a c , &) = (0.48, 0.1), D=0.2, 
i?=0.25, and £?bi q =0.025, where the energy unit is meV. 
This parameter set gives the frc-plane spin spiral with a 
wave number gf,=0.37r at low temperatures, which resem- 
bles the spin structure in TbMn0 3 (q 6 =0.297r) 

We trace dynamics of the Mn spins by numerically 
solving the Landau-Lifshitz-Gilbert equation using the 
fourth-order Runge-Kutta method. We derive an effec- 
tive magnetic field Hf s acting on the spin Si from the 
Hamiltonian % as Hf s = —dW/dSj. Considering the 



observed reduced Mn moment [2l|, we set the norm of 
the spin vector 1^1=1.4. The system used for calcula- 
tions is 40x40x6 in size with the periodic boundary con- 



we consider — E ■ p^ with 
This coupling effectively 
modulates the in- plane FM exchanges from J a bSi ■ Sj to 
(J a b — E ■ TVij)Si ■ Sj. Consequently, the applied E\\±a 
modulates the spin exchanges as shown in Fig.[ljb). Here 
|7r° | is calculated to be 3.5xl0~ 26 \xCm from the lattice 
parameters [22| and the observed ferroelectric polariza- 
tion P(~5000^C/to 2 ) for i?Mn03 with an up-up-down- 
down spin order [2j|. This means that E a =l MV/cm 
induces the modulation |A J a f,| = |£^ a 7r° |=0.022 meV. 

The infrared-absorption spectrum is calculated as the 
response to a weak <5-function pulse. (For technical detail, 
see Ref. [it})- The calculated spectrum is displayed in 
Fig. Die), which has two peaks at w=0.94 THz and 2.1 
THz, and reproduces well the experimental spectrum of 
TbMnOa. The Gilbert-damping coefficient is chosen 
to be 0.1 so as to reproduce the observed peak width, 
which guarantees the under-damped spin oscillations. 

Now we theoretically demonstrate switching of the spin 
chirality by the optical pulse. There are two kinds of 
6c-plane spirals with different spin chiralities, i.e., CW 
and CCW ones. Their spin chiralities C point in the 
—a and +a directions (C||— a and C||+a) so that they 
are referred to as 6c_ and bc + , respectively. Here the 
chirality C is defined as a sum of the local contribu- 
tions Ci_i-\-± = Si x Si+x and Ci^+y = Si x Si+y as 
C = 27v E.^'^+i+Cv+j)/^- Note that C hi+i is zero 
because of the AFM stacking in z-direction. The CW and 
CCW a6-plane spirals are also possible although they are 
slightly higher in energy than the 6c-plane ones without 
external fields. They have C||— c and C|j+c, and thus 
are referred to as ab- and a&+, respectively. 

Starting with frc_ with C||— a, we apply an intense 
pulse of E=[E a (t), 0, 0] along the a axis, where 



E a (t) = — Eq sin tot exp[— 



(t - t f 
2a 2 



(4) 



with to=l psec. Here the frequency uj is fixed at 2.1 THz, 
which corresponds to the higher-energy electromagnon 
peak, while the full width of the half maximum for the 
Gaussian envelope, 2\/2hi2<7, is taken to be 0.5 psec 
[see inset of Fig. [2]. In Fig. [21 we display calculated 
time evolutions of the a-, b- and c-axis components of 
C . When E$=+IA MV/cm, we find a reversal of C from 
be- (C Q <0) to bc+ (C a >0) as shown in Fig^a). This 
reversal takes place via ab + with C c >0 and C a ^0. In 
addition, as shown in Fig. EJb), when we apply a slightly 
weaker pulse of Eq=+13 MV/cm, a chirality flop occurs 
from be- (C||— a) to ab+ (C\\+c). These switchings oc- 
cur very fast, typically within 5-7 psec. 

To understand these phenomena, we first consider en- 
ergies of the four chirality states, i.e., ab± and bc±. When 
,E=0, the ab + and ab- are degenerate, and are higher in 
energy than the ground-state 6c-plane spirals. Applica- 
tion of i£||±a lifts this degeneracy, and the ab + (ab-) 
becomes the lowest in energy under the strong £^||+a 
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FIG. 2: Time evolutions of the spin chirality C=(C a , Cb, C c ) 
after application of the pulse with (a) _Eo=+14 MV/cm and 
(b) Eo=+13 MV/cm. Fig. 2(a) shows the chirality-reversal 
from be- {C a <0) to bc+ (C a >0), while Fig.2(b) shows the 
chirality flop from be- (C\\ — a) to ab+ (C\\+c). Insets show 
spin states before or after applying the pulse, and time profile 
of the applied pulse E a (t). 



(a) ab+ (CCW) 



(b) ab- (CW) 
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FIG. 3: Modulations of spin turn angles (thin red arrows) and 
the local spin chiralities ±Adj under .E||+a for (a) ab+ and 
(b) ab- states, (c) Arrangement of the c-axis components of 
DM vectors ±0^. 



(E\\—a). This can be understood as follows. As shown 
in Figs. (SJa) and GUJb), the modified in-plane FM ex- 
changes, J a b ± AJ ao , under E\\±a cause changes in the 
spin turn angles (thin red arrows) and hence staggered 
modulations of the local spin chiralities as Cij ± ACij 
(thick brown arrows). Under .E||-|-a, the modulations 
±ACij in the ab+ (ab-) are always antiparallel (paral- 
lel) to the staggered c-axis components of the DM vectors 
±df • [blue arrows in Fig. EJc)] . Since the DM coupling 
favors the antiparallel configurations of ±dj • and zhCi.j, 
the energy decreases (increases) in the ab+ (ab-). Note 




FIG. 4: (a) Schematic figure for time evolution of the po- 
tential structure in the 6 space during the chirality-reversal 
process for Eo=+14 MV/cm. Here 9 is the angle between spin 
chirality C and the a axis. For a time profile of the applied 
pulse with -Eo>0, see inset of (b). (b) Relationships between 
the chirality-switching processes and the sign of Eq. 



that this mechanism is distinct from the direct coupling 
between E and Pij given in Eq. (1), which is much weaker 
in the THz-frequency regime. On the other hand, the 
energies of bc± are not affected by E because the a-axis 
components of the DM vectors are alternately stacked, 
by which the DM energy always cancels out. 

In Fig. HJa), we illustrate schematic time evolution 
of the potential as a function of 9 during the chirality- 
reversal process. For time profile of the applied pulse 
with Eq>0, see the inset of Fig. 0|b). Here 8 is the an- 
gle between the chirality C and the a axis. At t=t\, 
the system is located in the minimum at 0=180° (be-). 
When E a (t)>0 as at t~t2, 0—90° (ab + ) becomes a new 
energy minimum, so that the chirality C starts rotating 
or the angle 9 starts decreasing towards this minimum. 
Importantly the chirality does not stop its rotation at 
9=90° immediately, but passes through that minimum 
or oscillates around it because of the inertial force. The 
DM interaction and the single-ion anisotropy originat- 
ing from the spin-orbit interaction make the rotation of 
spin-spiral plane massive, resulting in its inertial motion. 
Then the E a (t) becomes negative as at t=t% : which makes 
ab+ (0=90°) the highest in energy Consequently the sys- 
tem starts falling into the minimum at 9—0° (bc + ). At 
last (t=ti) : the chirality reversal is completed being set- 
tled in bc+ (6=0°). 

On the other hand, the chirality flop for Eq=+13 
MV/cm is a rather subtle process. If the system is within 
the domain of metastability of 8=90° when E a (t) be- 
comes almost zero at t=t^ the system can be trapped in 
the local minimum of a&+. Since the ab± are mctastablc, 
the system should decay into the ground-state bc + or bc- 
eventually due to thermal fluctuations. 

This ultrafast chirality switching is distinct from the 
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FIG. 5: (a)-(f) Color maps of calculated angles 9ij between 
the local spin chiralities and the a axis (see inset), which 
show the real-time dynamics of Cij. 



conventional polarization reversal, and is accompanied 
by the dynamical spatial-pattern formation, i.e., dynam- 
ical stripes of chirality domains, through the nonlinear 
photoexcitation of electromagnons. In Figs. H{a)-(f), we 
show the calculated real-time dynamics of the local chi- 
ralities Ci.j on the Mn-0 plane for the chirality-reversal 
process. Starting from be- (blue) at t=0, the optical 
pulse generates ab + (green) domains as at £=1.2 psec. 
Subsequently stripes of bc + (red) and be- (blue) domains 
emerge after the pulse ends at t=2 psec, and the chiral- 
ities Ci_j oscillate in each domain between 0^=0° and 
9ij=180° as seen in 2<i(psec)<3.4. Among the spirally 
ordered spins, those directing (nearly) along the propa- 
gation vector cannot flip and become nodes of the spin 
oscillations to form the chirality domains. At last (t=6 
psec) the system gradually gets settled in the 6c+ state. 

Finally we discuss conditions for the chirality switch- 
ing. First, the switching occurs only at a frequency of 
the higher-energy electromagnon resonance (a;~2.1 THz 
in the present case), and does not occur at the lower- 
energy peak (w~l THz). Second, we need a rather large 
peak height of the pulse, lEbl^lO MV/cm, in the present 
simulation. Experimentally maximum peak height ex- 
ceeding 100 MV/cm is available for 10-72 THz [2J], but 
below 3 THz, it reaches only ~1 MV/cm at present. We 
expect that the optical pulse with I-EqI^IO MV/cm at 
~2 THz will be realized in the near future. Besides, 
the threshold value of \E \ can be reduced if we prop- 
erly choose the target materials: For example, solid so- 
lutions Tbi_ x Gd T Mn03 locating near the boundary be- 
tween the ab- and 6c-plane spiral phases [25| are promis- 
ing candidates. Importantly there are optimal ranges of 
the electric strength \Eq\, and a larger \Eo\ cannot nec- 
essarily induce the switching. This can be understood 
as follows. To achieve the chirality reversal, for instance, 
from 6»=180° (6c_) to 9=0° (bc + ), the chirality vector 
oscillating around the energy minimum at 0=90° (ab+) 
should be in the range 0°<#<90° when E is reversed 
from E a (t)>0 to E a (t)<0 in order to fall into the an- 
other minimum at 6>=0° (bc + ) instead of 0=180° (6c_). 



This means that we need to adjust depth of the energy 
minimum of ab + at t=ti by tuning the strength of Eq in 
order to synchronize the timing between the chirality os- 
cillation and the reversal of E. Therefore, the switching 
processes show highly nonlinear behaviors with respect 
to strength and shape of the pulse. In addition, if we 
adopt a negative E in Eq. (4), the lowest- lying state at 
t=t2 becomes ab- with 0=270°. Then the chirality rever- 
sal occurs via ab- instead of ab+. The chirality flops to 
ab- from bc± become also possible for a slightly weaker 
\Eq\. Relationships between the switching processes and 
the sign of Eq are summarized in Fig. @|b). 

To summarize, we have theoretically studied the ul- 
trafast optical switching of spin chirality by exciting the 
electromagnons in the multifcrroic Mn perovskitc. We 
have revealed that the oscillating E component of the 
light activates the collective rotations of the spin-spiral 
planes with a THz frequency via the ME coupling, and 
their inertial motions result in chirality reversal or flop. 
By tuning strength, shape and length of the pulse, the 
spin chirality is shown to be controlled at will. 
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